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Abstract 



^ \ We introduce here the q-Laplace transform as a new weapon in Tsallis' 



arsenal, discussing its main properties and analyzing some examples. 
The q-Gaussian instance receives special consideration. Also, we de- 
rive the q-partition function from the q-Laplace transform. 
KEYWORDS: q-Laplace transform, tempered ultradistributions, complex- 
plane generalization, one-to-one character. 
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1 Introduction 



1.1 The Laplace transform 

The Laplace transform, introduced by Laplace in his work on probability 
theory, is a widely used integral transform with many applications in physics 
and engineering. It is a linear operator acting on a function f{t) with a real 
argument t) that transforms it to a function F{s) with complex argument s, 



This transformation is essentially bijective for the majority of practical uses; 
the respective pairs f{t), L{s) are matched in tables. The Laplace transform 
has the useful property that many relationships and operations over the 
originals f{t) correspond to simpler relationships and operations over the 
images L{s). The Laplace transform is related to the Fourier transform (FT), 
but whereas the FT transform expresses a function or signal as a series of 
modes of vibration (frequencies), the Laplace one resolves a function into 
its moments. Customarily, in speaking of the Laplace transform without 
qualification one means the unilateral or one-sided transform. The Laplace 
transform can be alternatively defined as the bilateral, or two-sided one. 




oo 



die **/(i); s complex. 
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by extending the limits of integration to be the entire real axis. If that is 
done the common unilateral transform simply becomes a special case of the 
bilateral transform where the definition of the function being transformed is 
multiplied by the Heaviside step function. 

1.2 q-Statistical theory 

Nonextensive statistical mechanics (NEXT) |3l HI |5], an extension of the 
standard Boltzmann-Gibbs (BG) one, is used in variegated scientific areas. 
NEXT is based on a nonadditive (although extensive [6]) information mea- 
sure characterized by the real index q (with q = 1 recovering the BG en- 
tropy). It has been employed in diverse scenarios such as cold atoms in 
dissipative optical lattices dusty plasmas [8], trapped ions [9], spin- 
glasses [in], turbulence [TT], self-organized criticality [12], high-energy ex- 
periments at LHC/CMS/CERN [I3| and RHIC/PHENIX/Brookhaven [H], 
low- dimensional dissipative maps [15], finance [TB], galaxies [TT], Fokker- 
Planck equation's applications [18j, etc. A typical NEXT feature is that can 
it can be advantageously cast using appropriate q-generalizations of standard 
mathematical concepts [i9\. Included are, for instance, the logarithm and 
exponential functions, addition and multiplication, etc. 
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1.3 Our aims 

Here we add, by recourse to ultradistributions (see Appendix) the Laplace 
transform tool to such armory. It is to be pointed out that, quite recently, an 
alternative form of the q-Laplace transform (qLP) has been advanced by Won 
Sang Chung [T], who uses for that purpose so-called q-sums, q-differences, 
q-products, and q-ratios, which renders the treatment rather abstract. 
Our qLP is, instead, based on the ordinary version of the four elementary 
arithmetic operations. More importantly, Wong Sang Chung's definition 
does not follow the tenets of the pioneer paper by Umarov-Tsallis-Steinberg 
[2] that introduced the qFourier transform. In particular, the function to 
be Wong-transformed does not contain a q-exponential argument, a crucial 
NEXT-aspect that must be respected so as to maintain the theory's consis- 
tency. 

The reader is advised to peruse the Appendix before embarking into our 
discussion below. 
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2 Laplace transform from Fourier one 

Before dealing with the q- Laplace transform it is convenient to show that 
the ordinary bilateral Laplace transform can be obtained from the complex 
Fourier transform. Thus, let (see Appendix for details and references) 

• be the space of distributions of exponential type, 

• U or Ui the space of tempered ultradistributions, 

• J- the Fourier transform connecting them, and 

• H the Heaviside step function. 
We have 

J- : Aoo ^ W, (2.2) 
reading {'^{k) is the imaginary part of k and ^{k) its real one) 

oo 

T{k) = H[^{k)] J f{x)e'^'' dx - H[-^{k)] J f{x)e"'^ dx. (2.3) 

-oo 

The associated inversion formula is 

fix)^^fnk)e-''^dk, (2.4) 

where the contour Fp surrounds all singularities of T{k) and runs parallel 
to the real axis from — oo to oo above it, and from oo to — oo below. Then, 



making the change of variables ik = —p, the bilateral Laplace transform is 
obtained (see [20] ) 

C-.A^^Ui (2.5) 

and given by 

oo 

£(p) = H[^{p)] J /(x)e-P^ dx - H[-^{p)] J f{x)e-P'' dx. (2.6) 

-oo 

Let us insist: Ui is the space of tempered ultradistributions and we have 
made the change of variables p = —ik. Now, the corresponding inversion 
formula is: 

f{x) = ^ {c{p)eP'dp, (2.7) 



27ri 

where the contour surrounds all singularities of C{p) and runs parallel to 
the imaginary axis from —ioo to ioo to the right of it, and from ioo to —ioo 
to the left. If we consider distributions of exponential type f{x) such that 
f{x) = for X < 0, we obtain the single Laplace transform 

oo 

C{p) = Hmp)] I f{x)e-P^ dx, (2.8) 







whose inversion formula is 



fix) = ^ I C{p)eP^ dp, (2.9) 



where is the right hand side of the path Tl- 



3 q-Laplace transform 



Let Q stand the space of functions of the real variable x that are parame- 
terized by a real parameters q. We have defined in (j2T] - [23]) the q- Fourier 
transform 



F -.Q — >U 



(3.10) 



as 



F(/)(fc, q) = Fik, q) = [H{q - 1) - H{q - 2)] x 

oo 

H[^{k)] [ /(x){l + i(l-g)A;x[/(x)](''-^)}^, dx- 



H[-'^{k)] j f{x){l + i{l-q)kx[f{x)]'^'''^^}^^ dx\ . (3.11) 



and its inverse transform as 



27r 



lim / F{k, q)5{q — 1 — e) dq 



e"^'^^ dk. 



Ff L 1 

Let Qi be the space of functions of the real variable x 



(3.12) 



nj = {f{x)/fix)eninnj}, 



(3.13) 



where 



= {/(x)//(x){l - (1 - g)px[/(x)](^-i)}T^ G C'[R+]; 



fix) > 0; < \x\'g{x)e^'';s,ae R+;pe Z; > 



and 



1 < g < 2} 



^7 = \f{x)/f{x){l - (1 - G C'[R-]; 



(3.14) 



fix) > 0; < \x\'g{x)e'''';s,ae R+;pe Z; 3?(p) < 



1 < g < 2} 



(3.15) 



Making again the change ik = —p we immediately obtain the bilateral q- 
Laplace transform L 

L-.rii — >Ui, (3.16) 



as 



L{f){p, q) = Lip, q) = [H{q - 1) - H{q - 2)] x 

oo 

H[^ip)] J /(x){l-(l-g)px[/(x)](«-i)}T^, dx- 



H[-^ip)] j /(x){l-(l-g)px[/(x)](<?-i)}i^dx 

— oo J 

The corresponding inversion formula is easily found from (13.121) 

r ^ 

fi^) = TT- f 1™ / Hp, Q)^i<l -l-t)dq 

l-Kl J e^0+ J 



L 1 



dk. 



(3.17) 



(3.18) 



If we consider / G f2/ such that f{x) = for x < we obtain the unilateral 
q-Laplace transform 

L{p,q) = [Hiq-l)-Hiq-2)]x 

oo 

H[^{p)] j f{x){l - (1 - g)pa;[/(x)](^-i)}T^ dx, (3.19) 



and its inversion formula 



2m 



lim / Lip, q)5(q — 1 — e) dq 

e-S>0+ 

1 



dk. 



(3.20) 



We consider now functions G i7/ depending on the parameter q with 
1 < g < 2. We can define the singular q-Laplace transform 



Lr : — > Ui, 



(3.21) 



as 



LR{f'„){p,q) = LR{p,q) = lim L(/ /)(p, g) = L{f >){p,q) 



q=q 



(3.22) 



As is the case of the Ft q- Fourier transform (see [23]), Lr is NOT one 
to one. To deal with such an issue we consider the set A// , given by 



and 



A,^, = [g, G n:/LR{g,){k) = LR{f.{k)] 



(3.23) 



(3.24) 



Introducing the equivalence relation 



9q' e A//^, 



(3.25) 



and the q-Laplace transform between equivalence classes 




(3.26) 



as 



LpR{Ajf^,){p,q') = LpR{p,q) = LR{f^>){p,q), 



(3.27) 



one finds that LpR is one to one between equivalence classes and is the analog 
of FjjTs (the Umarov-Tsallis-Steinberg q- Fourier transform [2j) for the one 
to one q-Fourier transform. 

4 Examples 

We illustrate here with some examples the preceding developments. As a 
first one we consider the q-Laplace transform of the Heaviside's step function 
f{x) = H{x). We have 



L{p, q) = [H{q - 1) - H{q - 2)]H[^{p)] / [1 + (g - l)px]^^]dx. (4.1) 



oo 







Suitably manipulating (14. ip leads to 



L{p,q) = [H{q-l)~H{q-2)] 



mm 

(2 - q)p 



(4.2) 
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For f{x) — H{—x) we have 



L{p, q) = [H{q - 1) - H{q - 2)]H[-^{p)] J [I + {q - l)px]^^]dx, (4.3) 

— oo 

and, as a result, 

L{p, q) = [H{q - 1) - H{q - 2)]^^^!^. (4.4) 

Taking now into account that H{x) + H{—x) = 1 we get for the q-Laplace 
transform of f{x) = 1 

L{p, q) = [H{q - 1) - H{q - (4-5) 

We evaluate now the transform of f{x) = q H{x), i.e., 

L{p, q, q) = [H{q - 1) - H{q - 2)]i^ x 

00 

J q'[l + {q- l)pxq'^'^-^'^]^^]dx. (4.6) 


One finds 

L(p, q, q) = [H{q - 1) - H{q - 2)f-^^^^^. (4.7) 

2. — q p 

If we consider now the Laplace transform of the previous function we face 
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(4.9) 



As a last example we evaluate the transform of the function 

{^Y ; X e[a,b] ; < a < b ; A > 

; X outside [a, b] . 
One has 

L{p, q) = [H{q - 1) - H{q - 2)]H[^{p)] x 

b 

J x-~^{l - (1 - g)pA^(^^^)x^-^(''~^)}^. dx (4.10) 

a 

Following the steps of a similar calculation made in [21] allow us to obtain 

Lip, q) = [H{q -l)-H{q- 2)]H\^{p)] x 



H{q-l)-H 



g- 1 + 



/3 



X 



(9-l)A^ 



■X 



g-2 / 1 

ai-^F 



{2-q)[{q-l)p\P] — 
2-q 1 



m - q) 

g-l'(g-l)[l-/?(g-l)]'g-l ' 1 - /3(g - 1) 



+ 



1 



„ I 1 



(1 - g)j9A'3(9-i)ai-/3(9-i) 
2-q 1 



+ 



m - q) 



g - 1 ' (g - 1)[1 - /3(g - 1)] ' g - 1 l-(5{q-iy 
1 



(1 - g)pA/^(9-l)6l-/3(3-l) 



if 



1 



-H{q-2) 
/3-1 



+ 

A^ 
(5q-2 



g-l'/3(g-l)-l'/3(g-l)-l' 
12 



(1 - ?)pA^(^-i)ai-^(«-i)) 



(l-g^A^Mfei-'^^^-^))}}. (4.11) 
and, if we take /3 = l/(g — 1), we get for out transform 

Lpr{p, q) = H[^{p)] [H{q - 1) - H{q - 2)] [1 - (1 - ?)pA]T^ . (4.12) 

5 Series expansion of the q-Laplace transform 

Consider the function 

{l-{l-q)px[f{x)r'}^^. 

Using the series expansions of the logarithm and the exponential function, 
we can write 

{1 - (1 - g)px[/(x)]^-^}^ = eT^Ml-(l-g)pa.[/(a.r-i} ^ 

1 - (^i}!!±i(g_l).(p^)n[_^(^)]n(,-l) 

e "=i = 



_L. J2 (q-l)r^|^px)ngn(q-l)lnf(x) 

e ' "=i 



g n=l m=0 



E E (px)"ln-[/(a=)](g-l)"+— 1 

.n=l m=0 
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E 



E ■ 

m=0 



-(pi;)"-'"+lln'"[/(i;)] 



-1)"} 



Let g{x,p,n) be given by 



m=0 



m! 



Then, 



or 



{1 - (1 - q)px[mr'} = = e™5o^(-^-)(^-^)" 



y~] g{x,p,n){q—l)" 

h{x,p,q) = e~^^e"=i 



Minding the expansion of the exponential function we have 



\n=l 



and, as a consequence. 



n=l 



where 



^ oo s— n+1 s— si— n+2 s— si— S2 Sn-2 — 1 

= 5] E E 

s=n si=l S2=l s„-i=l 

g(x, p, si)g{x, p, S2) • • • g{x, p, x 



g{x,p, S-S1-S2 Sn-i){q - ly 
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One can write the q-Laplace transform as 

L{p,q) = [H{q-l)-H{q-2)]x 

oo "1 

H[^{p)] J f{x)h{x,p,q) dx - H[-^{p)] J f{x)h{x,p,q) dx > . (5.8) 

-oo j 

6 The q-Laplace transform of the q-Gaussian 

Our purpose is to calculate the q-Laplace transform of the q-Gaussian. As 
this becomes a too complex task in the general case, we content ourselves 
with a first-order expansion in powers of g — 1. Accordingly, 

h{x,p, q) = e-P^[l + g{x,p, l){q - 1)], (6.1) 



with 



f 59x)^ 

g{x,p,l) ^ — pxln[f{x)]. (6.2) 



Then, up to first order we have for the q-Laplace transform 

L{p,q)^[H{q-l)-H{q-2)]x 

oo 

HMp)] / {l + -pxln[/(x)]| (q - 1)} /(x)e-^^ dx- 



\ 

H[-m] 1 1 1 + I ^ - ln[/(x)] }(?-!)} dx i . (6.3) 
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Let G{k) and G{k, (3) be given by 



G{p) = { H[^{p)] J /(x)e~P^ dx - H[-?ii{p)] I f{x)e-P'' dx 





G{p,P)= {Hmp)] j[f{x)fe-^'^dx-H[-^{p)] I [f{x)fe 





We write the q-Laplace transform in the form 



L{p,q) = [Hiq-l)-Hiq-2)]x 



G{p) + 



p2 Q2 



2 dpi^ip) + p--Gip, P) 



d_d_ 

dp 8(5 



-!)• 



Let f{x) be the q-Gaussian 



f{x) = CAl + {q -l)ax'Y^ 



where 



B 



1 1 1 

2' q'-l 2 



q 7^1, 



Using results of p5] we obtain 



G{p,q)=H[^{p)]C^ 2 



{q - l)ap^ 



2-9 
1-9' 
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and 



Vr ^ [ l-q' 



g' 9 



P+i-q 1 



H 



1-9 



^ 2 



N 



- l)ap_ 
P 



^"U(9'-i)« 



i7[-3?(p)]C"', 



p ; 13+1-q 



g' 9 ^ 

^ [(g'-l)«]T^ 



,8 + l-g 1 



H 



2 



N 



(1 - q')ap_ 
p 



X 



(6.11) 



w ^ v(i-^')«y v(i 

where H and N are the Struve and Neumann functions, respectively. The 
q-Laplace transform of the q-Gaussian is now 



L{p,q,q) = [H{q-l)-H{q~2)]x 



G{p,q') + 



^\G{p,q') + k^4;-G{p,q',P) 



(6.12) 



13=1 



If instead of using the bilateral q-Laplace transform we use the unilateral 
one, we should replace fl6.10p and (16. lip , respectively, by 

2 



G{p,q) = H[^{p)]C^ 



[(g'-l)a]^ 
17 



137 L(g' - l)ap. 



2-9 1 



X 



and 



13 \/'^ 



1 2 



H 



P 



-N, 



^ + 1 



{q - l)ap 

P 



1-9 



7 2 



X 



(6.14) 



7 The q-Laplace transform of the q- Gaussian 
for fixed q 

In this section we deal with the q-Laplace transform of the q- Gaussian for 
fixed q. We have 



Lpr{p, q) = H[^{p)] j C,[l + {q- l)ax']—^x 



|l + (1 - q)ikx jcjl + {q- l)^^^]^^} {q - l} dx- 



H[-^{p)] J Cg[l + {q-l)ax^]^x 

— oo 

1 

{l + (g-l)px{cjl + (g-l)Q;x']i^}(g- l}"' dx, (7.1) 
1 < q < 2. Simplifying terms we obtain 

oo 

Lpr(p, q) = H[?ft(p)] j C,[(q - l)ax' + (q - l)Cf^px + l]i^ dx- 



u 

H[-^{p)] J Cg[iq-l)ax^ + e-^{q~l)CfH:x + l]^dx. 



(7.2) 



Effecting tlie cliange of variables ^y{q'^-T)a x = y, the q-Laplace transform 
becomes 



a 



1 

1-9 



dy- 



H[-^{p)] 



+ 1 



1-q 



dy. (7.3) 



Using results given in [26j we see that (P(^ is the associated Legendre func- 
tion) 



2^r(i -2^)(^2 _ 



r(i-/x)r(-/x-z/)r(i/-/x + i) 



X 



(7.4) 



and we can write 



oo 

y (1 + 2te + t^)'^-^ t-^-''-^^ dt = T{~^i)2-^-\z^ - l)^P't^_,{z), (7.5) 



Let 7, /i be given by: 



7 



1 1 



/i 



1-q 2' 



(7.6) 



so that 



LpR{p,q) = C,^L4=2-^-\jy - l)f X 

v(g- 1)« 
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(7.7) 



which is the q-Laplace Transform of the q-Gaussian (on the complex plane) 
for fixed q. If instead of using the bilateral q-Laplace transform we use the 
unilateral one, we just obtain 

LpRip, q) = C,-^td=2~'^-~W - l)^//[5R(p)]P^_,(7P). (7.8) 
V(g - l)a 

8 The q-partition function 

As it is well known, the partition function of a system is the unilateral Laplace 
transform of the density of states [22]. Thus, the q-partition function should 
be defined as the q-Laplace transform of the density of states f{u), in which 
the complex variable is now B (in place of p) . We have 



If the density of states depends on q we can define it following the LpR 
definition given in section 2. 



Z{B,q) = [H{q-l)-H{q-2)]x 



oo 




(8.1) 







ZpR{B,q) = [Hiq - 1) - Hiq - 2)]x 



oo 




(8.2) 
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For example, if the density of states is a q-exponential 



l)au] 1-9 , 



(8.3) 



where a > 0, the q-partition function is 



ZpR{B,q)^[H{q-l) 



H{q-2)] 



B + a ' 



(8.4) 



An important case is f{u) — constant — H{x). We face a non-degenerate 
energy spectrum (the one dimensional harmonic oscillator, for instance). We 
have then. 



Conclusions 

We have here developed the q-Laplace transform, thus incorporating it to 
the TsaUis' arsenal. We have studied its main properties and analyzed some 
instructive examples. The particularly important case of the q-Gaussian has 
been discussed in some detail. Finally, we have derived the q-partition func- 
tion from the q-Laplace transform. As is also the case with the q-Fourier 
transform, we realize that the q-Laplace transform is essentially a transfor- 
mation between equivalence classes. 



Z{B,q)^[H{q-l)-H{q-2)] 



H[^+{B)] 
(2 - q)B ■ 



(8.5) 
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9 Appendix: Tempered ultradistributions and 
distributions of exponential type 

Many readers will surely benefit from a brief summary of the main properties 
of distributions of exponential type and of tempered ultradistributions. 
Notations. The notations are almost textually taken from Ref. [50]. Let 
(res. C") be the real (resp. complex) n-dimensional space whose points 
are denoted by x = (xi, X2, a;„) (resp z = {zi, Z2, Zn)). We shall use the 
notations: 

(a) x + y = {xi + yi,X2 + y2,---,Xn + yn) ; ax = {axi,ax2, ...,axn) 
(h)x ^ means Xi ^ 0, X2 ^ 0, x„ ^ 

n 

{c)x ■y=Yl XjVj 
i=i 

n 

(d)| x\=J2\xj\ 

Let N"' be the set of n-tuples of natural numbers. If p G N"', then p = 
{pi,P2, ...,Pn), and Pj is a natural number, 1 ^ j ^ n. p + q stands for 
{pi + qi,p2 + q2, ■■■,Pn + qn) and p^ q means pi ^ qi,p2 ^ g2, ^ qn- 

n 

x^ entails We shall denote hj \ p \= J2 Pj call the 

differential operator dP^+P^+-+P^ /dxiP^Xi^^ ...dxn^^ 

For any natural k we define x'^ = x\x\...x^ and /dx^ = d""^ / dx\dx\...dx\ 
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The space H of test functions such that e^l^l|D''0(x)| is bounded for any p 
and g, being defined [see Ref. ([30j)] by means of the countably set of norms 



p = sup ePl'^l 

0<q<p, X 



DU(x) 



, p = 0,l,2,... (9.1) 



The space of continuous hnear functionals defined on H is the space Aoo of 
the distributions of the exponential type given by ( ref. PU] )■ 

[e'=N/(x)] (9.2) 

where k is an integer such that k ^ and f{x) is a bounded continuous 
function. In addition we have H <Z S <Z S' <Z Aqo, where «S is the Schwartz 
space of rapidly decreasing test functions (ref|3T]). 
The Fourier transform of a function (p G H is 

oo 

^(z) = — [ J{x) e^^-^ dx (9.3) 
2n J 

— oo 

According to ref. [30], 0(z) is entire analytic and rapidly decreasing on straight 
lines parallel to the real axis. We shall call H the set of all such functions. 

n = T{n} (9.4) 

The topology in "H is defined by the countable set of semi-norms: 

, = sup|zh0(z)|, (9.5) 

zGVk 
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where Vk = {z = {zi, Z2, Zn) G :| ImZj \^ k,l ^ j ^ n} 
The dual of H is the space tl of tempered uhradistributions [see Ref. ([30] )]■ 
In other words, a tempered ultradistribution is a continuous hnear functional 
defined on the space H of entire functions rapidly decreasing on straight lines 
parallel to the real axis. Moreover, we have ncS CS CU. 

hi can also be characterized in the following way [see Ref. ([30] )]: let 
be the space of all functions F{z) such that: 

F{z) is analytic for {z G : \Im{zi)\ > p, \Im{z2)\ > p, \Im{zn)\ > 

p}. 

S)- is bounded continuous in {2; G : \Im{zi)\ ^ p,\Im{z2)\ ^ 
p, \Im{zn)\ ^ p}, where p = 0,1, 2, ... depends on F{z). 
Let n be the set of all 2;-dependent pseudo-polynomials, z G C"^. Then U is 
the quotient space 

C)-u = A^/n 

By a pseudo-polynomial we understand a function of z of the form 

^jG{Zi, Zj-i, Zj^i, Zn) with G{zi, Zj-i, Zj+i, Zn) G 

Due to these properties it is possible to represent any ultradistribution as 
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[see Ref. ([30] )] 

F(0) =< F{z),(j)iz) >= i F{z)(t){z) dz (9.6) 



r 

r = Fi U r2 U ...r„, where the path Tj runs parallel to the real axis from 
— cxD to oo for Im{zj) > C, C > p and back from oo to — oo for Im{zj) < 
~C < ~P- (r surrounds all the singularities of F{z)). 

Eq. fl9.6p will be our fundamental representation for a tempered ultradistri- 
bution. Sometimes use will be made of the "Dirac formula" for ultradistri- 
butions [see Ref. ([20])] 

oo 

F(z) = i / — dt (9 7) 

^ ' {2m)^ J (ti - Z,){t2 - Z2)...{tn -Zn) ^ ' ' 

— oo 

where the "density" f{t) is such that 

oo 

F{z)(t){z) dz= [ f{t)(l){t) dt. (9.8) 



While F{z) is analytic on F, the density f{t) is in general singular, so that 
the r.h.s. of (19. 8 p should be interpreted in the sense of distribution theory. 
Another important property of the analytic representation is the fact that 
on r, F{z) is bounded by a power of z 



\F{z)\<C\z\P, (9.9) 



where C and p depend on F . 
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The representation fl9.6l) implies that the addition of a pseudo-polynomial 
P{z) to F{z) does not alter the ultradistribution: 

{F{z) + P{z)}(f){z) dz = ^ F{z)(j){z) dz + ^ P{z)(f){z) dz 
r r r 

However, 

P{z)(j){z) dz = 0. 



As P{z)(f){z) is entire analytic in some of the variables Zj (and rapidly de- 
creasing), we obtain: 

{F{z) + P{z)}4>{z) dz = j F{z)(P{z) dz. (9.10) 
r r 
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